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Abstract 

q ■ From an experimental point of view, clear signatures of multifragmentation 

O ■ have been detected by different experiments. On the other hand, from a theo- 

retical point of view, many different models, built on the basis of totally dif- 
"q ' ferent and often even contrasting assumptions, have been provided to explain 

them. In this contribution we show the capabilities and the shortcomings of 
one of this models, a QMD code developed by us and coupled to the nuclear 
de-excitation module taken from the multipurpose transport and interaction 
code FLUKA, in reproducing the multifragmentation observations recently re- 
^ ! ported by the INDRA collaboration for the reaction Nb + Mg at a 30 MeV/A 

i | projectile bombarding energy. As far as fragment production is concerned, 

O ■ we also briefly discuss the isoscaling technique by considering reactions char- 

acterized by a different isospin asymmetry, and we explain how the QMD + 
| FLUKA model can be applied to obtain information on the slope of isotopic 

yield ratios, which is crucially related to the symmetry energy of asymmetric 
nuclear matter. 

vo : 

t^j- ! 1 Features of multifragmentation 

When studying heavy-ion collisions at non-relativistic energies, multifragmentation can be observed 

for the most central ones, in a range of projectile-ion bombarding energies from tens MeV/A up to a 
| few hundreds MeV/A, depending of the properties of the nuclei under consideration. Many issues of 

this phenomenon are still under discussion, in particular concerning the stage at which it occurs in the 
. evolution of a reaction, e.g. if a nuclear system undergoing multifragmentation is or not equilibrated, 

how a simultaneous break-up in multiple fragments can occur, and if the multifragmentation is the result 
5_i ■ of a phase-transition. 

According to the currently most believable scenario, during the overlapping stage of heavy-ion 
collisions (typical time ~ 100 fm/c) matter can undergo compression, leading to large excitation ener- 
gies. As a consequence, the blob of nuclear matter starts to expand and can go on expanding down 
to sub-saturation densities (p ~ 0.1 - 0.3 po, where po is the normal nuclear matter density) and reach 
temperatures ~ 3 - 8 MeV, where it becomes unstable and breaks up into multiple fragments. These 
conditions are typical of a liquid-gas coexistence region [1,2]. 

As already mentioned, one of the open issues is if equilibration is reached in these reactions. A 
statistical description of multifragmentation is based on this assumption. A dynamical description of 
multifragmentation instead is not based on this assumption. Difficulties in coming to a non-controversial 
conclusion are largely due to the fact that most of the experimental data refer to patterns of particles 
detected just in the last stage of the reactions, involving channels feeded by the sequential decays which 
heavily affect and modify the primary fragment distribution. 

Multifragmentation can be distinguished from other decay channels on the basis of the excita- 
tion energy: a typical scenario of small excitation energies (E < 2 - 3 MeV/A) is characterized by the 
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formation of a compound-like system and by its evolution through binary sequential decays (evapora- 
tion/fission), whereas for high excitation energies (E > 3 MeV/A) multifragmentation in a finite volume 
and a simultaneous break-up into multiple fragments can occur. The excitation energy is indeed related 
to the mass asymmetry (A pro j — A targe t): in case of symmetric central reactions the compression is 
responsible of the high excitation energy, whereas in case of asymmetric reactions only a partial com- 
pression can occur and a large part of the excitation energy appears in the form of thermal energy [3]. In 
all cases, multifragmentation is tipically driven by the following expansion. 

Due to its mentioned features, multifragmentation, occurring during the phase of expansion of the 
nuclear system formed by an ion-ion (central) collision, allows to study the nuclear Equation of State 
(EoS) at subnormal nuclear densities. In particular, it is possible to infer useful information concerning 
the symmetry energy and its density dependence, by investigating the isotopic yield distributions of the 
emitted fragments. The isoscaling technique, based on the analysis of isotopic yield ratios obtained in 
reactions with a different isospin asimmetry, has been developed with this purpose. 

After an overview of the models to study multifragmentation in Section and a brief presentation 
of the one used in this work in Subsection 12.11 examples of its application in the isoscaling technique 
and in the reconstruction of multifragmenting sources at energies of a few tens MeV/A are provided in 
S ubsection 13 . 1 1 and 13.21 respectively. Finally, our perspectives on further applications of our model are 
drawn in Section 0] 

2 Models to study multifragmentation 

One can distinguish between 

- Dynamical Models: some of them are 1-body approaches, inspired to the BUU/BNV/Landau- 
Vlasov transport theory. Alternatively, n-body approaches have been developed, such as the 
QMD/AMD/FMD. n-body approaches are very powerful in the description of the simultaneous 
break-up of a nuclear system in multiple fragments, since they preserve correlations among nucle- 
ons. 

- Statistical Models: they assume to work with an equilibrated excited source at freeze-out (thermal 
equilibrium). Taking into account that the nuclear system undergoes an expansion, leading to 
decreasing densities, down to subnormal values, the freeze-out [4] occurs when the mutual nuclear 
interaction among fragments can be neglected. Statistical models have been worked out both in 
the grand-canonical framework (see e.g. Ref. [5]) and in the micro-canonical framework. The 
most widespread among the last ones is the SMM [1,6,7] and its modifications ISMM [8] and 
SMM-TF [9]. 

We emphasize that the onset of multifragmentation according to dynamical models is different 
from the description of multifragmentation according to statistical models. In fact, in the statistical mod- 
els a source in thermal equilibrium is assumed to fragment. This means that memory effects concerning 
how the source has been originated are neglected. On the other hand, in the dynamical models multifrag- 
mentation is a fast process: the involved nucleons have not the time to come to equilibrium. Fragments 
originate from the density fluctuations (nucleon-nucleon correlations) due to collisions in the ion-ion 
overlapping stage, which survive the expansion phase (memory effects). The chemical composition of 
hot fragments is expected to play a role in helping to disentangle the nature (dynamical / statistical) of 
the multifragmentation mechanism [10]. 

Different models reproduce different features of the collisions with different success. A mixed 
model, inspired to the QMD dynamical approach to describe the fast stage of ion-ion collisions and to 
a statistical approach to describe the further decay of the multiple primary excited fragments produced 
by QMD down to their ground state, has been used to obtain the results presented in this work. Due to 
the crucial role of dynamics, as supported by our results, in the following we mainly concentrate on the 
description of the dynamical aspects of multifragmentation. 
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2.1 QMD/AMD approaches 

In these microscopic models a nucleus is considered a set of mutually interacting nucleons. The prop- 
agation of each nucleon occurs according to a classical Hamiltionian with quantum effects [11]. In 
particular, nucleons are described by gaussian wave packets. Each of them moves under the effects of 
a potential given by the sum of the contribution of all other nucleons (2-body effects). Furthermore, 
when two nucleons come very close to each other, they can undergo elastic collisions (nucleon-nucleon 
stochastic scattering cross-sections) with Pauli blocking. 

A proper treatment of antisymmetrization is implemented in AMD [12]. On the other hand, QMDs 
do not provide any antisymmetrization of the nuclear wave-function. An approximate effect can be ob- 
tained through the inclusion in the Hamiltonian of a Pauli potential term, or through the implementation 
of specific constraints. 

Still open questions in molecular dynamics approaches concern the functional form of the nucleon- 
nucleon potential (each working group who developed a molecular dynamics code has its preferred 
choice of terms), the potential parameters and their relation to the nuclear matter EoS. Nowaday, many 
groups prefer parameter sets leading to a soft EoS. Anyway, there are open questions concerning the 
symmetry term [13, 14]. In particular, a stiff dependence for this term means that the symmetry energy 
always increases with increasing densities. On the other hand, a soft dependence means that the sym- 
metry energy decreases at high densities. At present, a stiff dependence seems more reliable than a soft 
one. Many uncertainties come from the fact that our observations are mainly based on symmetric nu- 
clear matter (N / Z ~ 1) near normal nuclear density (~ 0.16 fm~ 3 ), since it is difficult to obtain highly 
asymmetric nuclear matter in terrestrial laboratories. On the other hand these studies are crucial to under- 
stand features of astrophysical objects (such as neutron star formation and structure), where conditions 
of extreme neutron-proton asymmetry can be present. 

Other open issues concern the gaussian width, the use of in-medium nucleon-nucleon cross- 
sections instead of free nucleon-nucleon cross-sections (in QMD the free choice is usually implemented, 
whereas in the AMD the in-medium choice has been implemented), the question of how long the dynam- 
ical simulation has to be carried over and the problem of the development of a fully relativistic approach 
(on the last point see e.g. Ref. [15]). 

A QMD code has been developed by us [16] in fortran 90. It includes a 3-body repulsive potential 
and a surface term (attractive at long distances and repulsive at short distances). Pauli blocking is imple- 
mented by means of the CoMD constraint [17]. Neutron and proton are fully distinguished by means of a 
simmetry term and an isospin dependent nucleon-nucleon stochastic scattering cross-section. The kine- 
matics is relativistic and attention is paid to the conservation of key quantities (total energy/momentum, 
etc.) in each ion-ion collision. Simulations are performed by means of our code from the ion-ion over- 
lapping stage up to t ~ 200 - 300 fm/c (fast stage of the reaction). The description of the de-excitation 
of the excited fragments present at the end of the fast stage is obtained through the coupling of our 
QMD with the statistical model taken from the PEANUT module available in the FLUKA Monte Carlo 
code [18-21] in a version for the g95 compiler. Up to now, the QMD + FLUKA interface has been tested 
in the collisions of ions with charge up to Z=86 (radon isotopes), providing interesting results (see e.g. 
Ref. [22] and references therein). 

3 Results 

3.1 Isospin dependence in fragment production: application of the isoscaling technique 

The isoscaling technique, already mentioned in Section [Q is based on ratio of yields taken in multifrag- 
mentation reactions with similar total size, but different isospin asymmetries (N — Z) I (N + Z) [23,24]: 

R 21 (N, Z) = Y 2 (N, Z)/Y X {N, Z) = Const exp(A coeff N + B coeff Z) . (1) 
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The numerator of this formula refers to the yield of a given fragment (N, Z) obtained from a neutron rich 
nucleus-nucleus reaction system, whereas the denominator refers to the yield of the same fragment from 
a neutron poor (more symmetric) reaction at the same energy. A coe f f is related to the symmetry energy 
and is increasingly larger for couple of reactions with increasingly different isospin composition N/Z. 

In particular, we have considered the neutron rich systems Ar + Fe (N/Z = 1.18) and Ar + Ni 
(N/Z = 1.13) with respect to the neutron poor system Ca + Ni (N/Z = 1.04). Among other authors, 
these systems have been previously studied by [25] (see also Ref. [26,27]). 

Isotopic yield ratios for light fragment (Z < 8) emission have been obtained from our QMD + 
FLUKA simulations for the couple of reactions Ar + Ni / Ca + Ni and Ar + Fe / Ca + Ni at 45 MeV/A 
projectile bombarding energy. When plotted in the logarithmic plane, isotopic yield ratios for each fixed 
Z turn out to be approximately linear, with a slope given by a A coe ff, as expected from Eq. (Q]). As for 
the isoscaling parameter A coe ff, our simulations give the following insights: 

- The results of our analysis are quite sensitive to the number of isotopes included in the linear fit, at 
fixed Z (i.e. to the goodness of the gaussian approximation to the fragment isotopic distribution). 

- A coe f j differs with Z, in agreement with [23], which claims that isoscaling is observed for a variety 
of reaction mechanisms, from multifragmentation to evaporation to deep inelastic scattering, with 
different slopes in the logarithmic plane. 

- A^ff is larger for the couple of reactions with larger difference in the isospin compositions 
Qfx/Zx - N 2 /Z 2 ). 

- Our average values A coe f f = 0. 18 for Ar + Ni / Ca + Ni and A coe f f = 0.31 for Ar + Fe / Ca + Ni 
are larger than the experimental values [25], but the comparison is not so meaningful, since it is 
largely affected by the fact that we include fragments emitted in all directions in our preliminary 
analysis, whereas in the experiment only fragments emitted at 44° were selected. 

- A coe f f turns out to be affected by the choice of the impact parameter and decreases significantly 
when selecting only the most central events. 

- A^f f,hot at tne en d °f the overlapping stage can be larger than A coe ff at the end of the full 
simulation by no more than 20%, at least for the reaction systems under study. 

As far as the emissions at preequilibrium are concerned, our simulations lead to the following 
results: 

- For central collisions of Ca + Ni the yield of emitted protons turns out to be larger than the yield 
of emitted neutrons by 20%. For central collisions of Ar + Ni and Ar + Fe, on the other hand, the 
yield of emitted protons turns out to be lower than the yield of emitted neutrons by 10 - 15%. 

- For each of the three systems under study, the fragment asymmetry of the liquid phase (Z/A)u q at 
the end of the preequilibrium stage turns out to be lower than the corresponding value at t = 0, in 
qualitative agreement with the AMD simulations [25]. 

- No traces of isospin fractionation appear, expected indeed for systems with an higher N/Z content 
(e.g. 60 Ca + 60 Ca). 

The dependence of our results on the projectile bombarding energy is currently under study, by 
considering the same reactions at different bombarding energies. 

3.2 Multifragmenting source reconstruction in Nb + Mg reactions at 30 MeV/A 

Multifragmentation has been observed in Nb + Mg reactions at a 30 MeV/A projectile bombarding en- 
ergy in an experiment performed at the INDRA detector by the INDRA + CHIMERA collaborations [28]. 
Event selection has been performed, according to experimental cuts on the momentum along the beam 
axis, p z ^det > 0.6p 2; tot, and on the angular acceptance of the INDRA detector, 4° < 9 < 176°. The 
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selected events have then been assigned to different regions, corresponding to portions of the plane iden- 
tified by the total transverse energy and the total multiplicity of charged particles detected in each event. 
Three regions have been singled out this way, as shown in Fig. 2 of Ref. [28]. We have applied the 
same selection procedure by implementing proper cuts and filters on the simulated events obtained by 
our QMD + FLUKA. The selected theoretical events are plotted in Fig. [Q which can be directly com- 
pared with Fig. 2 of Ref. [28] and turns out to be in good agreement. The events plotted in the Tl region 
(red) are the less dissipative ones (more peripheral collisions), whereas the events in the T3 region (blue) 
correspond to more dissipative (central) collisions. 

Total detected charged particle multiplicity vs total detected transverse E 

30 | , , , , , 1 

25 
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1 1 1 1 1 1 1 

50 100 150 200 250 300 

E_trasv_det (MeV) 

Fig. 1: Multifragmentation of Nb + Mg at 30 MeV/A: event selection and identification of different regions Tl 
(red), T2 (green) and T3 (blue) by our QMD + FLUKA simulations. Each point corresponds to a different ion-ion 
reaction event in the plane identified by the total multiplicity of detected charged particles and the detected total 
transverse energy. 



For each of the three regions, average values of interesting quantities have been obtained both 
in the experiment and in the theoretical simulations. Our results, concerning the transverse energy, the 
multiplicity of charged particles, the velocity and the charge of the biggest residual averaged over all 
events belonging respectively to the Tl, T2 and T3 regions are shown in Table [TJ As far as the average 
transverse energy and multiplicity of charged particle are concerned, the results of our simulations turn 
out to be in good agreement with the experimental data, within the experimental uncertainties, in the 
region Tl and T2, corresponding respectively to peripheral and semiperipheral collisions, whereas in 
case of central collisions the theoretical average transverse energy underestimates the experimental one 
and the theoretical multiplicity of charged particles slightly overestimates the experimental result. On 
the other hand, as for the properties of the largest residual, the results of the theoretical simulations show 
good agreement with experimental data especially for the most central collisions, belonging to the T3 
region, whereas for the more peripheral ones the theory overestimates the velocity and the charge of the 
largest residual. These results, considered all together, seem to point out to the fact that in the experiment 
the interacting nuclei are slightly more stopped than in the simulation. 

Furthermore, the considered experiment aim at the recostruction of the properties of the so-called 
source, the blob of matter formed by compression in the ion-ion overlapping stage, which undergoes 
multifragmentation. Since the experiment detects final cold fragments, i.e. fragments in their ground 
state after the de-excitation, a procedure has been established to reconstruct the properties of the source 
by using the observed properties of the final fragments and emitted protons. In particular, the source is 
isolated by a selection in parallel velocity of different fragments (velocity cuts), by considering different 
cuts in different regions. The velocity cuts implemented are summarized in Table 2 of Ref. [28]. As 
far as the proton are concerned, the parallel velocity cut is fixed to 3 cm/ns, whereas for increasingly 
heavier fragments the velocity cuts are fixed to increasing values. The velocities of the emitted fragments 



Table 1: Multifragmentation in Nb + Mg at 30 MeV/A: average values of interesting quantities in the regions Tl, 
T2 and T3. 



Region < E trasv > < M tot > 


Vmaxres -- > 


^ Z-maxres 


Experiment (from Ref. [28]): 


Tl 72.9(10.1) 8.1(0.8) 


6.6 (0.2) 


34.1 (2.3) 


T2 120.9(10.6) 11.2(0.9) 


6.4 (0.2) 


31.9 (2.7) 


T3 176.4(13.8) 13.4(1.0) 


6.3 (0.2) 


30.5 (2.8) 


Theoretical simulations (QMD + FLUKA de-exc.): 


Tl 74.8 8.4 


7.0 


38.8 


T2 115.8 12.0 


6.6 


35.7 


T3 155.1 15.3 


6.2 


31.1 



are easily obtained even from our simulation by QMD + FLUKA, so it is possible to apply the same 
procedure for the reconstruction of the source properties even in case of our simulation. As an example, 
the velocities of the emitted protons obtained by our simulation for events in each of the three regions 
are shown in Fig. |2j by plotting their perpendicular component v perp vs. their component along the 
beam axis v par . This figure can be compared with Fig.3 of Ref. [28]. The vertical line in each panel 
corresponds to the v par cut implemented in the reconstruction of the source. 



region T1 region 12 region T3 




-5 5 10 15 20 -5 5 10 15 20 -5 5 10 15 20 

v z ( cm/ns ) v z ( cm/ns ) v z ( cm/ns ) 

Fig. 2: Multifragmentation of Nb + Mg at 30 MeV/A: v perp vs. v par for protons emitted in the region Tl (left 
panel), T2 (central panel) and T3 (right panel) , respectively, as obtained by our QMD + FLUKA simulations. 
Each point in each panel correspond to a different emitted proton. The vertical lines correspond to the velocity cuts 
implemented for the reconstruction of the multifragmenting sources. 

Since the experiment is able to detect the charge of the emitted fragments but not their mass, the 
velocity cuts can be directly used just to obtain the charge of the source Z s . To calculate the mass of the 
source A s a further hypothesis is needed. The author of Ref. [28] assume that the source has the same 
isotopic ratio as the projectile, i.e. A s /Z s = A pro j /Z pro j. Source properties in the three regions Tl, 
T2 and T3, as reconstructed both from the experiment and from our simulation, are shown in Table [2] 
Since the theoretical model allows to simulate even the process of source formation in a straightforward 
way, the properties of the source can be directly obtained before its de-excitation and break-up into 
multiple fragments, without using an a-posteriori reconstruction based on velocity cuts. If we identify 
the source with the biggest fragment present just at the end of the QMD simulation, we obtain a very 
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good agreement with the experimental results of Ref. [28], especially in the region Tl and T2, even if 
the experimental results are based on the a-posteriori reconstruction, as can be inferred from Table |2] 
On the other hand, if we use a reconstruction procedure analogous to the one adopted by the authors of 
Ref. [28], we overestimate the size of the source, especially for the most central collisions. Finally, as 
for the average multiplicity of IMF fragments (Z > 3) subsequently emitted from the source, we obtain 
good agreement with the experiment in all regions. 



Table 2: Multifragmentation in Nb + Mg @ 30 MeV/A: reconstruction of the source properties in the regions Tl, 
T2 and T3 



Experiment (from Ref. [28]): 


Region 


<z s > 


< A s > < M p > 


<M a > 


< M frag > 


Tl 


40.7 (2.0) 


91.2(4.7) 2.0(0.6) 


1.1 (0.5) 


1.2 (0.4) 


T2 


42.8 (2.1) 


96.0 (4.9) 2.7 (0.7) 


1.8 (0.6) 


1.4 (0.3) 


T3 


45.1 (2) 


101.3(4.6) 3.1(0.7) 


2.5 (0.7) 


1.6 (0.3) 


Theoretical simulation: 


Region 


<z s > 


< A s > < M p > 


< M a > 


< M fra g > 


Source properties reconstruction from final secondary fragments (QMD + FLUKA de-exc): 


Tl 


42.6 


96.2 2.2 


0.4 


1.0 


T2 


45.1 


101.5 3.5 


1.6 


1.25 


T3 


48.5 


109.5 4.2 


3.8 


1.5 


Source properties at the end of the QMD simulation (primary fragments): 


Tl 


41.0 


91.0 






T2 


43.3 


96.6 






T3 


47.5 


106.1 







4 Conclusions and perspectives 

The QMD model developed in Milano and coupled to the de-excitation module of the Monte Carlo 
FLUKA code has been used to study reactions between ions of intermediate mass which exhibit multi- 
fragmentation features. The results presented in this paper are encouraging, and can be further refined 
by more precisely investigating up to which extent the statistical de-excitation process from FLUKA 
modifies the pattern of primary fragments originated dynamically by QMD, and how the results of the 
simulation change when the time of the uansition from the dynamical description of the nuclear system 
to a statistical description is modified. 

Further studies at non relativistic energies that we are going to perform with our theoretical simu- 
lation tool concern: 

- the isospin distillation effect: it occurs in the multifragmentation of charge-asymmetric systems, 
and leads to IMF fragments (liquid) more symmetric with respect to the initial matter, and light 
fragments (gas) more neutron rich. This effect is related to the density dependence of the symmetry 
energy. 

- The bimodality in the probability distribution of the largest fragment as a function of the mass 
number A max of the largest fragment, as a signature of a phase transition. Experimental data on 
this effect have been obtained by the CHIMERA collaboration (see e.g. Ref. [29]). 

- (Complete) fusion cross-sections (this kind of analysis has already been performed by other groups, 
e.g. by means of the ImQMD model [30]). 
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